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A hierarchical ordering is demonstrated for the periodic orbits in a strongly coupled multidimen- 
sional Hamiltonian system, namely the hydrogen atom in crossed electric and magnetic fields. It 
mirrors the hierarchy of broken resonant tori and thereby allows one to characterize the periodic 
orbits by a set of winding numbers. With this knowledge, we construct the action variables as 
functions of the frequency ratios and carry out a semiclassical torus quantization. The semiclassical 
energy levels thus obtained agree well with exact quantum calculations. 
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I. INTRODUCTION 

The hydrogen atom in crossed electric and magnetic 
fields is among the paradigmatic examples of strongly 
coupled multidimensional systems. During the last two 
decades, a large number of experimental and theoretical 
investigations (eg., 0, 0, IE 0, IE IE S H, H, E3, E3, E3, Q ) 
have been devoted to the intricate physics of this seem- 
ingly elementary nonintegrable system. In addition to its 
inherent interest, the hydrogen atom in crossed fields can 
be used to study gener ic p henomena such as Arnold diffu- 
sion 0, monodromy or Ericson fluctuations 
Its deeper understanding also provides a solid basis to 
explore the behavior of confined electrons in condensed 
matter physics, such as in excitons |14| and quantum dots 
|l5j . Yet, in spite of the large amount of work spent upon 
it, the overall phase space structure of the crossed-fields 
hydrogen atom still defies a complete understanding. 

One of the most prominent features of a dynamical 
system are its periodic orbits (POs). While individual 
POs can yield deep insights into the local dynamics, in 
their entirety they provide an appropriate tool to under- 
stand the global structure of a multidimensional phase 
space |l6j . Their knowledge is crucial in many appli- 
cations of classical dynamical systems such as astron- 
omy 0, 0, , particle accelerators flp a nd fluid dy- 
namics, e.g. statistics of turbulent flow [2(|. The semi- 
classical quantization of the classical structures is an in- 
valuable tool for the description of molecular vibrations 
pH 12^ . l23| , chemical reactions [24| or the spectra of Ry- 
dberg atoms 

In this work we establish an organizing principle for 
POs in the crossed-fields hydrogen atom. It allows the 
identification of winding numbers for every PO that 
originates from the breakup of a torus. Conversely, 
through this classification the POs elucidate the higher- 



dimensional structures of the phase space and allow one 
to characterize their topology. They enable us to recon- 
struct the hierarchy of broken tori with sufficient pre- 
cision to carry out an Einstein-Brillouin-Kcllcr (EBK) 
torus quantization pi IM I2I l2ll29l| . 

The electron motion of a hydrogen atom exposed to 
an electric field F in the x direction and a magnetic field 
B in the z direction, is governed by the Hamiltonian, in 
atomic units l30l . 



TT 1 2 1 B , 

H= 2 P — r +^( p y x 



B 2 

Px y) + —(x 2 +y 2 )-Fx. (1) 



Here r = (x, y, z) 

P = {Px,Py,P. s 



y/x 2 + y 2 + z 2 . 



are the usual Cartesian coordinates, 
the conjugate momenta, and r = 
By virtue of the scaling properties j^H 



of the Hamiltonian Q , if all classical quantities are mul- 
tiplied by suitable powers of the magnetic field strength, 
the dynamics can be shown not to depend on the energy 
E and the field strengths B and F separately, but only 
on the scaled energy E = B~ 2 / 3 E and the scaled electric 
field strength F = B~ 4 / 3 F. In particular, coordinates 
scale according to f = B 2 l z v and the classical actions 
obey S = B X / 3 S. In this paper, we present results for a 
scaled electric field strength F — 0.5 and two scaled ener- 
gies E = —1.5 and E = —1.4 slightly below and slightly 

above the classical ionization threshold E\ = —2\[f, re- 
spectively. 

The organization of this paper is as follows. In Sec. [H] 
we give a brief general review of the phase space struc- 
tures in Hamiltonian systems that are relevant to our 
investigations. Sec. IIIII explains the dynamical princi- 
ples underlying the organization of periodic orbits and 
summarizes the organization of periodic orbits of the hy- 
drogen atom in crossed electric and magnetic fields. In 
Sec. HVI we describe this organization in detail and show 
how it can be exploited to assign winding numbers to the 
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periodic orbits. Sec. El presents the calculation of action 
variables from the information thus obtained. A semi- 
classical Einstein-Brillouin-Keller quantization is carried 
out in Sec. IVII and conclusions are given in Sec. IVIII 



II. PHASE SPACE STRUCTURES IN 
HAMILTONIAN SYSTEMS 

Among all Hamiltonian systems, those that exhibit the 
most regular dynamics are the integrable systems, where 
almost every bounded trajectory is confined to an in- 
variant torus. At the other extreme, in ergodic systems 
almost every trajectory comes arbitrarily close to every 
energetically allowed point in phase space. In this work 
we treat the crossed-fields hydrogen atom in a parameter 
range where the phase space structure is dominated by 
the remnants of invariant tori. Nevertheless, the external 
fields are chosen too strong for perturbation theory to be 
a reliable tool, and we will consider the full nonintegrable 
Hamiltonian Q). 

By definition, a Hamiltonian system with / degrees of 
freedom is integrable if it possesses / constants of motion 
in involution. The regular level sets of these constants 
are /-dimensional tori if they are compact [3^, |3^| . For 
singular values of the constants, tori of lower dimensions 
can arise. Tori of different dimensions form a hierarchy 
in which the tori of dimension N are organized around 
the tori of dimension N — 1. Such hierarchies have been 
observed in a variety of different systems (see, e.g., 0, 
EH- EH)- I n the following sections, we will demonstrate 
how this hierarchy of tori in an integrable limit can give 
rise to a hierarchy of isolated POs in a nonintegrable 
system and how these POs can in turn be used to gain 
information on the original hierarchy of invariant tori. 

In an integrable system, action-angle variables (I, 8) 
with the following properties can be introduced. The 
angles 8 determine the position on an individual torus. 
They increase linearly with time 



8 = Ut + 8n 



(2) 



with a constant frequency vector u> and initial conditions 
#o- The conjugate action variables I are constants of mo- 
tion and characterize the invariant tori. They are given 
by 



1 

2^ 



ptfq, 



(3) 



where is the loop on the torus obtained as the angle 
0i varies from to 2ir with all other angles held fixed. 

The action variables 10 form the basis for a semi- 
classical torus quantization |25l 1 2(1 12 7L l28l |29| . Find- 
ing them explicitly, however, is a highly nontrivial task. 
Due to its significance in applications such as the cal- 
culation of molecular vibration spectra, the construction 
of the invariant tori has attracted considerable attention 



In an integrable system, the Hamiltonian can be writ- 
ten as a function of the action variables only, and the 
frequencies u) in @ are given by 

- w (4) 

If the ratios of any two of the u>i on a torus are rational, 
that torus carries periodic orbits and is called a resonant 
torus. In this case, a set w of integer winding numbers 
can be found so that 



0J\ : U>2 ■ 



Lu n — wi : u> 2 : . . . : w n . (5) 

Each winding number u>i specifies the number of rota- 
tions along the direction of the fundamental loop ji that 
the POs execute before repeating themselves. 

Angle coordinates 8 on a given torus can be defined 
in various ways. Apart from a choice of origin, which 
is inconsequential, any two angle coordinate systems are 
related by a linear transformation 

8' = M 8 (6) 

where M is an n x n integer matrix with unit determi- 
nant [3£|. The action coordinates must be transformed 
according to 



I' = (M ) ■ I 



(7) 



whereas the winding numbers on a rational torus trans- 
form as 

w' = M ■ w. (8) 

This freedom to choose different action-angle coordinate 
systems, together with the intimate connection between 
coordinate systems and winding numbers, will play a cru- 
cial role in our investigation of the phase space topology. 

In a non-integrable system, a resonant torus breaks 
up into isolated POs Ha H3, E3| . We use these POs as 
representatives of the torus they stem from and call them 
A^-torus POs, where N is the dimension of the original 
torus. According to Kolmogorov-Arnold-Moser (KAM) 
theory 32, 33], most nonresonant tori remain intact in 
a near-integrable system and break up only gradually. 
They are interspersed with the isolated POs in the same 
way as resonant and nonresonant tori are interspersed 
in an integrable system. Therefore, POs can be used to 
investigate the structure of the surviving tori. 

According to the Poincare-Birkhoff theorem jU H3, 
l39l |. a 2-torus in a two-degree-of-freedom system breaks 
up into an even number of isolated POs, usually two. 
One of the POs is stable ("elliptic"), the other is unstable 
("hyperbolic"). Because both POs stem from the same 
torus, their periods and actions differ only slightly. 

For systems with three degrees of freedom, Kook and 
Meiss [40j were able to derive from symmetry consider- 
ations that a 3-torus breaks up into four isolated POs 
that usually represent the four possible stability combi- 
nations: one PO elliptic in both degrees of freedom (ee), 
one PO hyperbolic in both degrees of freedom (hh) and 
two POs elliptic in one degree of freedom and hyperbolic 
in the second (eh and he). 
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Figure 1: Periodic orbits at E = —1.4 and F = 0.5. The 
FPOs S and their repetitions are shown with black circles, 
(a) the 2-torus POs T% (red plus symbols) and T% (green 
triangles), (b) the T% and their 3-torus partners (blue 
crosses), (c) the Tj 1 and their 3-torus partners T% (magenta 
diamonds) . 



III. ORGANIZATION OF PERIODIC ORBITS 

Our investigation of the POs in the crossed-fields hy- 
drogen atom starts with a numerical search for periodic 
orbits. This search, which is described in detail in Ap- 
pendix^ produces a long and unstructured list of POs. 
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Figure 2: The two planar FPOs S + (green) and S (blue) at 
E — —1.5 (a) and E = —1.4 (b), where S~ is already strongly 
deformed. The position of the nucleus is marked by a dot. 



Several such lists are available in the literature (see, e.g., 
El ) but no comprehensive ordering scheme for 
POs in three degrees of freedom has been proposed so far, 
and a priori it is not even clear if one exists. However, 
in Fig.^a clear structure underlying the family of POs 
becomes apparent: The periods and actions of most POs 
fall into well-separated series. (A few orbits that arise 
from secondary bifurcations and do not fit this pattern 
were omitted from the figure.) In a recent publication 
[4lj we argued that the POs arise from the breakup of 
invariant tori (i.e., they are A-torus POs in the termi- 
nology of Sec. [HJ and that the series structure provides 
evidence for a hierarchical ordering of POs that reflects 
the hierarchy of invariant tori in an integrable limit of 
the dynamics. In the followin g;, w e will briefly summa- 
rize the conclusions reached in |4l| , describe in detail the 
computational procedures that justify them and demon- 
strate how the ordering of POs can be used to gain in- 
sight into the topology of higher-dimensional phase space 
structures. In the course of the exposition, the nomen- 
clature used in Fig. ^ will be made clear. 

The three shortest or fundamental periodic orbits 
(FPOs) of the crossed-fields hydrogen atom were found 
in jj, [5j . They do not arise from the breakup of a higher- 
dimensional torus and therefore represent the 1-torus 
POs in the hierarchy. Two of these, labeled S + and S~ 
in p| , serve as organizing centers for the entire hierarchy. 
Both of them are planar. They are depicted in Fig. For 
very low energies, S + and S~ are shaped like Keplerian 
ellipses pj. At the energy levels chosen in this work, S + 
maintains this form, whereas S~ is strongly deformed as 
shown in Fig. [21 In addition, above the ionization thresh- 
old the FPO S~ , although still stable, is surrounded by 
a phase space region of ionizing trajectories P. IH l42| . 

The mechanism that enables the FPOs to organize the 
hierarchy of A-torus POs around them can best be illus- 
trated if the discussion is temporarily restricted to the 5- 
y plane perpendicular to the magnetic field. This plane 
forms an invariant two-dimensional subsystem. Its dy- 
namics is therefore accessible to an investigation using 
a Poincare surface of section such as Fig. 01 The plot 
shows the two FPOs S + and S~ , each at the center of an 
island of stability. They are surrounded by POs of larger 
periods. In an obvious way, the latter can be thought of 
as arising from the breakup of resonant invariant tori and 
being interleaved with surviving KAM tori as described 
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Figure 3: Periodic orbits in the Poincare surface of section 
y = for the planar subsystem, E — —1.5, F = 0.5. Non- 
fundamental POs are labeled by their winding ratios. 



in the previous section. Therefore, the POs can be la- 
beled with winding numbers and action variables that 
stem from the original family of invariant tori, and thus 
provide access also to the distribution of the remaining 
invariant tori, which are much harder to compute. 

An integrable Hamiltonian system with / > 2 degrees 
of freedom possesses an entire hierarchy of invariant tori 
of dimensions up to /. In every dimension the fully res- 
onant tori break up into isolated POs. A generic higher- 
dimensional Hamiltonian system will therefore exhibit a 
hierarchy of A-torus POs, where at each level the family 
of A-torus POs is organized by the underlying family of 
(A — l)-torus POs. This scenario is entirely analogous 
to the simple two-dimensional example of Fig. [21 but be- 
yond two degrees of freedom we cannot rely on Poincare 
surface of section plots to diagnose the situation. 

In Fig. all series end at lines that are formed by 
the FPOs and their repetitions. This observation gives 
us a first hint that the FPOs indeed serve as organizing 
centers for the families of longer POs even in the full, 
three-dimensional dynamics. This hint, however, is weak, 
and it is desirable to characterize the relation between the 
FPOs and the families they organize in a more detailed 
way. 

In Ref. ^3 ' we proposed three quantitative criteria to 
this end: (i) The stability angles <j>f of the FPOs S ± 
(i.e. the phase angles of the unimodular eigenvalues of 
their stability matrices) describe the rotation that each 
FPO imposes upon its neighborhood. The winding ratios 
of the 2-torus POs converge toward 4>f /2ir as the FPOs 
are approached, (ii) In the same limit, the action variable 
corresponding to the degree of freedom along the FPO 
converges to the action of the FPO, and (hi) the action 
variable for the motion transverse to the FPO, which is 
given by the area the (original) invariant torus encloses 
in the Poincar plane, tends to zero. We will use these 
criteria, which are derived in a two-dimensional setting, 




Figure 4: The hierarchy of iV-torus POs in the crossed-fields 
hydrogen atom at E — —1.5, F = 0.5 (and at E = —1.4, 
F — 0.5, except for the collapse on S~). The two families of 
3-torus POs T£ and T3 1 collapse onto their 2-torus partners 
Tf and T£ as ^3 — > 0. The T2 in turn collapse onto the 
FPOs S + and S~ . For a correct description of the collapse 
the appropriate coordinate system must be used: I\ — » for 
S~ and J( — > for S + . The coordinate change is described 
by the topological invariants M p and M n . 

more generally to characterize the relationship between 
the A-torus POs at different levels of the hierarchy. 

The entire hierarchy of POs that will be reconstructed 
with the help of these criteria is illustrated in Fig.0J We 
distinguish between two families T p and T n of A-torus 
POs that appear as series with positive (p) or negative 
(n) slope, respectively, in Fig. ^ Each of these contains 
a family T 2 p,n of 2-torus POs and a family T£' a of 3-torus 
POs. Since POs which are remnants of the same torus 
have nearly identical scaled actions S and periods T, each 
point in Fig. ^ represents one specific broken torus. The 
FPOs S + and S^serve as organizing centers for both 
families T 2 P and T 2 n of 2-torus POs. The POs located in 
the x-y symmetry plane (and shown in Fig. |3J) form the 
family T p . The 2-torus POs themselves are found to be 
limiting cases of two families T p and Xjf of 3-torus POs. 

To justify this description of the hierarchy, it is nec- 
essary to assign winding numbers to individual POs and 
calculate the corresponding action variables, so that the 
three criteria can be applied. These calculations will be 
presented in Sec. IIVI and [V] It will also be shown that 
the freedom of choice in the definition of action-angle co- 
ordinate systems that was explained in Sec. ITT1 and that 
at first sight seems to impede the application of the crite- 
ria can in fact be used to assign winding numbers to each 
torus in the coordinate system that is best adapted to the 
local dynamics. The way how different local coordinate 
systems transform into one another then contains global 
information about the topology of the families of invari- 
ant tori that is represented by the topological invariants 
MP and M n in Fig. [H 

IV. ASSIGNMENT OF WINDING NUMBERS 

A. The 2-torus POs T 2 P and T 2 n 

Winding numbers for the 2-torus POs T p and T-f can 
be assigned as follows: counting the series in Fig. ^ (a) 
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yields the first winding number w\. Counting the POs 
within one series from bottom to top yields the second 
winding number W2- The first series must be assigned 
i«i = 1 because the repetitions of the orbits in the left- 
most series occur in all higher series and the fc-fold repeti- 
tion of a PO results in another PO with winding numbers 
multiplied by k compared to the primitive PO. 

Using the same multiplication principle, we can per- 
form a first consistency check on this numbering scheme 
because it allows us to derive a pattern for the values 
of W2 across different scries. All POs in the first series 
are obviously primitive, whereas in higher series primi- 
tive and repeated POs are interleaved. As illustrated in 
Fig. [5] the winding numbers of all orbits in the higher 
series can be deduced from those of the first series. The 
numerical results shown in Fig. confirm this pattern. 
Only the base value w® to be assigned to the PO with 
the lowest action in the first series remains undetermined 
from these considerations. It will turn out that w® can- 
not be determined uniquely. Instead, different subsets of 
POs suggest different values of the winding number W2- 
As will be explained below, this non-uniqueness reflects 
the freedom of choice of different action-angle coordinate 
systems and allows one to assign winding numbers to 
each torus in the coordinate system that is best adapted 
to the local dynamics. 

To determine w® we use a Fourier series expansion of 
the time series of the Cartesian coordinates, with the 
frequencies u> measured in units of the fundamental fre- 
quency tDo = 2tt/T. For a PO with period T, these 
Fourier spectra do not show peaks at all multiples of the 



fundamental frequency, but only at those that correspond 
to the windin g nu mbers and some of their integer linear 
combinations [331 • For the planar Tf family it suffices 
to consider x(t), since the spectra for y(t) show peaks at 
identical positions. The first peak at u) = 1 • Qq in the 
example of Fig. \7\ corresponds to the series number w\ 
of the POs. It is found in the same position for any PO 
in a given series. The second major peak corresponds 
to the second winding number w 2 • Fig. |7fa) shows the 
spectrum for the PO with the lowest action in the first 
series, where this second peak is located at ui = 2 • Qq. 

We can therefore set w® = 2, which according to Fig.[S] 
determines the values of W2 for all orbits. However, as 
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Figure 5: The pattern of primitive and repeated POs for the 
first four series of 2-torus POs derived theoretically. The first 
series Wi = 1 contains only primitive POs. The second series 
wi — 2 shows an alternating sequence of primitive POs (odd 
W2) and repetitions of the POs in the first series (even W2). 
(It cannot be derived from the multiplication principle if there 
is an additional primitive PO at the lower or upper end of 
the second series.) In the third series there are blocks of 
two primitive POs (W2 indivisible by 3) separated by 3-fold 
repetitions of the first series POs (u>2 divisible by 3). The 
fourth series contains primitive POs for odd values of W2, 4- 
fold repetitions of the first series for positions with W2 divisible 
by 4 and 2-fold repetitions of the second series for W2 divisible 
by 2 but not by 4. The scheme can be extended to arbitrarily 
large values of Wi. 
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Figure 6: The pattern of primitive and repeated POs for the 
first four series of the Tj* (a) and the T£ (b). POs in the first 
series and their repetitions are shown with red circles, POs in 
the second series and their repetitions with green diamonds 
and POs from the third and fourth series in blue squares and 
black triangles, respectively. The sequence of primitive and 
repeated POs within the individual series should be compared 
to the pattern shown in Fig. Q>] which was derived from the 
multiplication principle. The numerical results agree perfectly 
with the theoretical predictions. 
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Figure 7: Fourier spectra (absolute values of Fourier coeffi- 
cients) of x(t) for some POs of the Tf family. Orbits are 
labeled by winding numbers calculated in two different coor- 
dinate systems (see text). The right column shows the orbits 
in configuration space and the position of the nucleus. 



Fig. E{b) and (c) show, the peak corresponding to this 
value of W2 is not always dominant and can indeed be 
very weak if W2 is large. Instead, a strong peak arises at 
a different position and suggests choosing a different set 
of winding numbers that has tu!; = 4 an d is connected to 
the previous system by the transformation 



w' 2 



Wi, 

2wi + w 2 - 



(9) 



Both choices are equally viable. They correspond to dif- 
ferent choices of angle coordinates on the original tori. 
As the Fourier spectra illustrate, the "unprimed" coor- 
dinate system is well adapted to the dynamics of orbits 
with low values of W2 , whereas the "primed" coordinates 
describe the dynamics of orbits with high w%. 

In view of Fig. (a) it is evident that POs with low 
W2 lie close to S~, while POs with high W2 are close to 
S + . The transition indicates that the dominant modes 
of the dynamics change from one end of the series to the 
other: close to S~ , these normal modes are the motion 
parallel and transverse to the FPO. As the distance from 
S~ grows, the coupling of the normal modes increases, 
and as S + is approached, the normal modes imposed by 
that FPO, which are different, gain dominance. In this 
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Figure 8: Fourier expansions of x(t) (red) and z(t) (green) for 
POs from the Tj 1 family. Peaks not associated with winding 
numbers are located at integer linear combinations thereof 



way, the two FPOs that bound the T p family impose two 
different systems of angle coordinates and 8' on the 
tori in their neighborhood. The transformation between 
the two coordinate systems is a topological invariant that 
characterizes the family of 2-torus POs. In our case, it is 
given by Eq. J!|J| and characterized by the matrix 



M p = 



1 

2 1 



(10) 



As anticipated in Sec. [nj M p is an integer matrix with 
unit determinant. 

Using a similar method, one can assign winding num- 
bers to the POs of the family. These POs are not 
restricted to the x-y symmetry plane, they are truly 
three-dimensional in Cartesian coordinate space, so that 
a Poincare surface of section plot similar to Fig. [3] can- 
not be obtained. Nevertheless, the Fourier spectra of x(t) 
and z(t) provide the information needed for a complete 
assignment of winding numbers, as shown in Fig. [S] It is 
again found that the dominant normal modes are differ- 
ent for POs at the lower end of the series, close to S~ , and 
at the upper end of the series, close to S + . The invariant 
AI n = Al p that characterizes the transition between the 
two limits takes the same value for both families. 



7 



B. The 3-torus POs If and T 3 n 

Apart from the 2-torus POs, Fig. ^ (b) and (c) show 
POs that are generated in the destruction of 3-tori. They 
are arranged in the same series as the 2-torus POs, and 
each 3-torus PO is intimately related to a 2-torus "part- 
ner" with almost identical period and action. Such part- 
ners have identical winding numbers w± and u>2, but 3- 
torus POs possess a third winding number u>3, which 
distinguishes different 3-torus PO partners of the same 
2-torus PO. It manifests itself in an additional peak in 
the Fourier spectra and can therefore be assigned by a 
straightforward extension of the technique used to clas- 
sify the 2-torus POs. For the T 3 n as in Fig. M (a) and 
(b) the situation is unambiguous since there is only one 
strong additional peak. The location of this peak pro- 
vides the winding number W3. For POs of the family 
(see Fig. |23(c) and (d)), however, two additional peaks of 
comparable magnitude arise, giving us two possibilities 
to assign the third winding number. In Fig. 0(d), the two 
possibilities are W3 — 5 and w' 3 = 1. As for the 2-torus 
POs, the two systems of winding numbers correspond to 
different angle coordinate systems on the original tori. In 
our case, the transformation between the two systems is 
given by: 

IU3 = wi - w 3 . (11) 

Note that there is no obvious relation between the 
classification by winding numbers that becomes possible 
though the present analysis and the geometrical appear- 
ance of the POs in configuration space. 



C. The stability of the POs 
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The stability of a PO in a Hamiltonian system with 
three degrees of freedom is characterized by a 4 x 4 sym- 
plectic stability matrix, which describes the linearized 
dynamics transverse to the PO [3^, 0, E3| • Because the 
stability matrix is real and symplectic, with each eigen- 
value A its inverse 1/A and its complex conjugate A* must 
also be eigenvalues. Therefore, if the eigenvalues of the 
stability matrix are different from ±1, they must belong 
to either (i) an elliptic pair e ±lip of complex conjugate 
eigenvalues with unit modulus, (ii) a hyperbolic pair A, 
1/A of real eigenvalues, or (iii) a loxodromic quartet A, 
1/A, A*, 1/A* of complex eigenvalues. A PO is stable 
if and only if all eigenvalues of its stability matrix have 
unit modulus. 

On an A-dimensional resonant torus in an integrable 
system, POs occur in continuous A-parameter families. 
They therefore have marginal stability, i.e., all four eigen- 
values of their stability matrix equal one. In a non- 
integrable system only a small number of isolated POs 
remains. The eigenvalues of their stability matrix occur 
in elliptic (e) or hyperbolic (h) pairs. In our numeri- 
cal studies we did not find any loxodromic quartets for 



Figure 9: Fourier spectra of x(t) (red) and z(t) (green). The 
second winding number is given in the unprimed coordinate 
system w° = 2. (a) A 2-torus PO in the T 2 n family and (b) 
its 3-torus partner from T£. (c) A 2-torus PO from the TJ 5 
family and (d) its 3-torus partner from T%. 



A-torus POs. In some cases our search algorithm finds 
a large number of POs all of which originate from the 
breakup of the same torus and have almost marginal sta- 
bility. In these cases, the splitting of the original torus 
is so small that we cannot resolve isolated POs within 
the given numerical precision. In the other cases, when- 
ever our numerical PO search allows us to identify all 
members in a quadruplet resulting from the breakup of 
a 3-torus, we find that it contains all four stability com- 
binations ee, eh, he and hh, in accordance with |4Cj . 

The 2-torus POs occur as doublets, of which Fig. ITUI 
shows an example. One pair of stability eigenvalues is 
approximately the same for both members of a doublet. 
It corresponds to the motion transverse to the original 
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Figure 10: Elliptic (green) and hyperbolic (blue) POs with 
winding numbers 1:2 (a) in the TJ 5 family and (b) in the Tj 1 
family. 

torus. With very few exceptions, the dynamics in this 
direction is stable. In the direction along the original 
torus one partner is elliptic while the other is hyperbolic 
in accordance with the situation of the Poincare-Birkhoff 
theorem |3£j in two degrees-of- freedom. Figs. [7||H1 El show 
elliptic (e) POs for the T 2 and completely elliptic (ee) for 
the T3 in those cases where isolated POs can be identified. 

D. Symmetry considerations 

The Hamiltonian possesses three discrete symme- 
tries, namely the reflection symmetry with respect to the 
x-y plane (z-parity), the y-parity with additional time re- 
versal and the combination of both. Due to the z-parity, 
the x-y plane constitutes a two-degree-of-freedom sub- 
system. Any PO outside this plane must either be itself 
invariant under the symmetry transformation or possess 
a partner related to it via the symmetry transformation. 

We find the following connection between winding 
numbers and symmetry: A PO in the T$ family is sym- 
metric under the z-parity transformation if and only if 
its winding numbers W\ and w 2 are even. In the T 2 and 
T3 1 families, the POs with w\ odd and w 2 even are sym- 
metric under z-parity, while all others are not. For both 
families of 3-torus POs the third winding number W3 has 
no influence on the symmetry properties. 

V. CALCULATION OF ACTION VARIABLES 

In this section we compute the scaled action variables 
I for each family of POs. 

The total action of a PO can be expressed in terms of 
the individual actions I: 

S = wl(w). (12) 

For any multiple a of the PO, this becomes: 

aS = aw • I (aw) , (13) 

with the same action variables I (aw) = I(w)as before. 
Thus, the I for the T3 families can be written as functions 
of any two frequency ratios (which by virtue of Eq. (J5J 



equal the ratios of the winding numbers) . Here we choose 
lK,™ 2)U ,3)=l(— ,— V (14) 

\w 2 W 2 J 

To calculate the action variables I numerically we select 
three POs close to each other in the (wi/w 2 , W3/W2) 
plane, so that the action variables are approximately con- 
stant in the small triangular area in between. Under this 
assumption, we obtain a system of equations 

§m = w^h+wfh+wfh, (15) 

which can be solved for I\, i 2 , and I3. These values are 
then assigned to be the function values at the barycenter 
of the triangle. The process is repeated for all POs in the 
corresponding family and the results are used to calcu- 
late an interpolation function using a modified Shepard's 
method Q. 

For the 2-torus POs we follow a similar procedure, ex- 
cept that we need to consider only a single winding ratio 
W\/w2- At the energy E = —1.5, we thus obtain the 
results shown in Fig. ^2 (a), which also displays the ac- 
tions and stability angles 0i i2 of the FPOs. The limiting 
values for high frequency ratios of the T 2 and T 2 fami- 
lies coincide with 4>i/2ir and </> 2 /27r of S~ , respectively. 
At the same time, the action variable I 2 converges to- 
ward the action of the FPO S~ , whereas the action vari- 
able It vanishes. According to the three criteria listed 
in Sec. IIIII we can thus conclude that S~ serves as an 
organizing center for both families T 2 and T 2 and that 
the action variable I2 and I± correspond to the normal 
modes along and transverse to S~ . (Notice that in three 
degrees of freedom a stable FPO has two pairs of unimod- 
ular stability eigenvalues, so that the collapse scenario of 
Fig. |3 can take place in two transverse degrees of free- 
dom independently, giving rise to two families of 2-torus 
POs.) For the planar family T 2 , this result merely con- 
firms the conclusion that we could already draw from the 
Poincare surface of section plot in Fig. By contrast, 
the family T 2 , although it consists of 2-torus POs, is not 
contained in any two-dimensional subsystem that could 
be described without an intimate study of the dynamics. 
Its analysis is therefore beyond the reach of a Poincare 
plot. Nevertheless. II If a) demonstrates that the relations 
of the families T 2 P and T 2 n to the FPO S~ are entirely 
analogous and that S~ organizes T 2 just as much as T 2 . 

As described in Sec. lIV Al the winding numbers w used 
in Fig. ^] (a) are ill-suited to describe the approach to 
S + , and the winding numbers w' given by © should be 
used instead. Indeed, none of the actions Jj and J 2 tends 
to zero in the limit of low frequency ratios. Nevertheless, 
S + can be located in Fig. ^2 (a) if its stability angles cj>' , 
which naturally arise in the w' system, are transformed 
to the w system just as if they were regular frequency 
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Figure 11: The actions 7i and 7~2 for the T% (red plus symbols) 
and T£ (green triangles), (a) Actions at E = —1.5, F — 0.5 in 
the unprimed action-angle coordinate system with w^- The 
winding numbers are suitable to describe the approach to S~ , 
where I\ tends to zero. The stability angles of the FPO S + 
have been transformed using Eq. 1171 . (b) The same data in 
the primed coordinate system with = 4. The situation 
around the FPOs is reversed, (c) Actions at E = —1.4, F — 
0.5 (unprimed coordinates). The T^'" do not collapse onto 
S~ due to the ionizing region around that FPO. (d) Actions 
at E = —1.4, F — 0.5 (primed coordinates). 



ratios by applying the inverse transformation of (j5J): 



f'2 



2w[ 



and therefore 



2tt 



(16) 



(17) 



The transformed stability angles <j> in Fig. ^2 ( a ) 
demonstrate that the T|"' n emanate from S + just as they 
emanate from S~ in the limit of high winding ratios. 
This result can be illustrated more clearly if the actions 
are calculated in the w' coordinate system, as shown in 
Fig. lllT b). Here, the situation around the FPOs reverses: 
J{ tends to zero as S + is approached and remains finite 
in the vicinity of S~. (Note that by virtue of Eq. Q 
there can exist coordinate systems in which the action 
variables take negative values.) In Fig. Illf b). the true 
stability angles of S + can be used to indicate the lower 
limits of the winding ratios, whereas the stability angles 
of S~ need to be transformed with the inverse of Eq. (|17fl . 

This symmetry demonstrates that locally the collapse 
of the Trf' 11 onto S + and S~ looks the same when de- 
scribed in suitable local coordinates. The nontrivial 
topology of the families Tj 5 and Tf that is represented 
by the invariants M p ' n becomes visible only if the en- 
tire families, including both limiting FPOs, are studied, 
cf.Fig.il 



Fig. ^3 (c) and (d) shows the situation at E — —1.4 
and F — 0.5, above the ionization saddle point. Here S~ 
is surrounded by an area filled with ionizing trajectories 
El H> This ionizing region prevents the T%' n from 
collapsing onto S~ . Yet, they still emanate from S + . 

Fig. El ( a ) shows the frequency ratios of 2- and 3-torus 
POs. For the 2-torus POs, the missing ratio CJ3/CJ2 is re- 
placed with the stability angle that describes the dynam- 
ics transverse to the original torus, normalized by 2ttw2- 
For the FPOs this angle is given by <fri — <p2- The sur- 
roundings of a long stable PO can rotate multiple times 
around the PO during one period. Stability angles, how- 
ever, can be calculated only modulo 27r. Therefore an 
integer number of 2-7T needs to be added to the stability 
angles of long POs. These stability angles arise in the 
u>3 coordinate system and thus need to be transformed 
with the inverse of Eq. i|ll[l . From Fig.^](a) it becomes 
apparent that the T3 approach the T2 in a suitable limit 
just as the T 2 approach the FPOs. 

Frequ ency maps such as Fig. I12f have been found in, 
e -g-; [lfll.l4fij to capture the essential dynamics of a multi- 
dimensional system. Their most prominent features are 
the resonance lines, each of which is given by a resonance 
condition of the form 



(18) 



with integer coefficients rrij. Their significance arises 
from the fact they they represent cantori that act as 
partial barriers to phase space transport. Diffusive tra- 
jectories such as those used in 0, ^EJ will stick to the 
cantori for a long time before leaving their neighborhood 
and thereby cause the prominence of the corresponding 
resonance lines in the frequency map. 

The periodic orbits shown in Fig. 1121 are the remnants 
of fully resonant tori, on which / — 1 independent res- 
onance conditions of the form Ijl8(l are satisfied. (They 
can be derived from the condition Q). 3-torus POs, 
therefore, will be located at the intersections of two res- 
onance lines in the frequency map. Since POs are non- 
wandering, they belong to the rigid background thr ough 
which the diffusion of generic trajectories takes place jlOj . 
As Fig.ll2lclearlv shows, the set of POs presents the web 
of dominant resonances just as clearly as the set of diffu- 
sive trajectories that is customarily used. 

The three action variables that characterize the 3-torus 
POs are displayed in Fig. ^] for the T$ family as func- 
tions of the two independent frequency ratios 6J1/C02 and 
They were calculated from Eq. 1)15(1 as described 
above. The action variables I\ and I2 obtained from the 
3-torus POs converge towards those found for the 2-torus 
POs, while -Z3 tends to zero as the lower boundary is ap- 
proached. We thus find all three criteria of Sec. [ffl] sat- 
isfied in this higher-dimensional situation: The tranverse 
frequency ratio 0J3/UJ2 of the tends to the transverse 
stability angle of the T^, the longitudinal action variables 
l\ and I2 of the 3-tori approach those of the 2-tori, and 




13 




0.1 0.15 0.2 0.25 0.3 0.35 



Figure 12: (a) Frequency map for JV-torus POs at E — —1.5, 
F = 0.5. The 3-torus POs T 3 P (blue crosses) and T3 1 (magenta 
diamonds) are bounded by the corresponding 2-torus POs Tlf 
(red plus symbols) and T 2 n (green triangles). The 2-torus POs 
in turn are bounded by the FPOs S ± (black dots), (b) The 
T% at E = —1.4, F — 0.5 with some of the most prominent 
resonance lines highlighted. They intersect in the points Pi = 

\T> l) ■ 



the transverse action variable I3 vanishes. We can there- 
fore conclude that the Tj 5 family serves as an organizing 
center for the T$, and we can identify the degree of free- 
dom corresponding to I3 as being transverse to the family 
of 2-tori. Similar results can be obtained for the collapse 
of the T§ family onto T<f, thus giving a full justification 
for all relations depicted in Fig. 

Because of these relations, additional data points at 
the lower boundary in Fig. 1131 can be obtained from the 
T% (and are included in the figure). For these, the missing 
frequency ratio is replaced with the stability angle 

according to Fig. ^1 the action variables 1\ and I2 are 
taken from Fig. 1111 and I3 is zero. 



Figure 13: Action variables for the T£ family at E = —1.4, 
F — 0.5. At the lower boundary, where the T£ are bounded 
by the T|", 73 tends to zero while 7i,2 tend toward the values 
obtained from the 2-torus POs. For clarity the dotted line /3 
from Fig. ^] (b) has been subtracted on the vertical axis. 



VI. TORUS QUANTIZATION 

The action variables calculated in the preceding para- 
graph provide the basis for an Einstein-Brillouin-Keller 
(EBK) torus quantization of the hydrogen atom in 
crossed electric and magnetic fields. The way of record- 
ing a quantum spectrum that is best suited to semiclassi- 
cal investigations of atomic spectra is scaled-energy spec- 
troscopy, which has for that reason been used in most 
experimental and theoretical work 0, 0, HH • It offers the 
advantage that the underlying classical dynamics does 
not change across the spectrum. A scaled spectrum con- 
sists of those values of the scaling parameter B~ x / 3 that 
characterize the quantum states for given scaled energy 
E and scaled electric field strength F. 

The phase space volume filled by 3-torus POs from the 
family is considerably larger than the volume filled by 
T3 1 , as can be seen by the number of respective POs for 
example in Fig.^ Thus, the largest part of the spectrum 
is also obtained from the quantization of the actions be- 
longing to T 3 P . 

The EBK quantization condition reads [M El IH |H 

Ii = j> pdq = 2irh (rn + , i = l, ...,/, (19) 

where the ji are the fundamental loops on the torus used 
in the definition @ of the action variables, rii are inte- 
ger quantum numbers and cii the corresponding Maslov 
indices. As shown above, one can obtain the action in- 
tegrals from PO data without having to determine the 
paths ji. In our case, we calculate the scaled action vari- 
ables U which need to be rescaled to obtain the true 
actions Ii = IiB^ 1 ^ 3 . Thus, the quantization condition 
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reads, with h = 1 in atomic units: 

I i (^)B-^ = 2.(n i + %). (20) 
\w 2 w 2 J V 4/ 

We derive the Maslov indices a,i from two simple con- 
siderations: In the limit of small I\ and I3, the action 
variable I 2 and its conjugate angle correspond to the mo- 
tion along the elliptic FPOs, which clearly has rotational 
character and requires a 2 — 0. The other two modes 
are transverse to this fundamental motion and are thus 
expected to have vibrational character, which leads to 
oei = 03 = 2. 

Inverting Eq. I|20|l to calculate B 1 / 3 gives: 

fl-i/a = + (21) 

1 y W2 ' W2 J 

A quantization condition for the scaled spectrum is ob- 
tained by observing that B -1 / 3 must take the same value 
for all three degrees of freedom. We can thus calculate 
the frequency ratios corresponding to the state with given 
quantum numbers n±, n 2 , 713 from the following set of 
equations 
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Table I: Semiclassical and exact eigenvalues of the scaling pa- 
rameter at E — —1.4 and F = 0.5. 

as in Table [I] yields the following correspondence between 
the torus quantum numbers n\, n 2 , n 3 and the quantum 
numbers n, q, k: 




and then compute the value of the scaling parameter 
from (|2T1) . 

Fig. [21 shows the EBK spectrum obtained from the 
quantization of the T|\ The low- lying states are labeled 
with the semiclassical quantum numbers (m, n 2 , 723). 
The manifolds of constant n 2 , corresponding to the quan- 
tization of I 2l determine the principal series discernible 
in the spectrum. The n2-manifolds overlap for n 2 > 4. 
For given n 2 , n\ ranges from to 1n 2 — 2, giving a total 
number of 2n 2 — 1 subseries. The number of states within 
one subseries is n 2 — Int [^1 . Subseries are identified by 
a constant value of m + TI3. The states within one sub- 
series are labeled by n 3 , beginning with for the highest 
state. 

Fig. El demonstrates that the spectrum derived from 
the quantization of the Tj 5 is in very good agreement 
with the exact quantum spectrum obtained as in, e.g., 
[T2I lit! |47| . The energy levels can be characterized in 
terms of three quantum numbers n, Qjk, which were 
first introduced in perturbation theory |48|: The princi- 
pal quantum number n identifies the principal series. The 
second quantum number q runs from — (n — 1), . . . , +(n — 
1) with increasing energy and corresponds to the sub- 
series. Finally, k counts the states within a subseries and 
runs from to n — \q\ — 1. A line- by-line comparison be- 
tween the semiclassical and the exact quantum spectrum 



n = n 2 , 

q = (m + n 3 + 1) - n 2 , (24) 
k = n 3 . 

Since I 2 corresponds to the fastest rotation, it comes as 
no surprise that the semiclassical n 2 represents the prin- 
cipal quantum number n. 

For large n 2 , the frequency ratios calculated from 
Eqs. (1221) and 12-il) for the lowest peaks in the central 
subseries (ni + n 3 rj n 2 ) lie outside the frequency range 
covered by the Tj\ These peaks belong to the TJ 1 family. 
However, the boundary of frequency ratios, especially in 
regions where few POs are available for the calculation 
of the action variables, is sometimes not clear-cut, which 
leads to a fuzzy boundary between the states of T£ and 
those of T£. Fig. H31 shows a spectrum for a single sub- 
series where the lowest states are taken from the quan- 
tization of the TJ 1 . The quantization of the T£ requires 
ot\ = and yields a correspondence between torus quan- 
tum numbers and traditional quantum numbers that is 
different from 1241) . The comparison with the exact quan- 
tum spectrum in Fig. 1 151 shows excellent agreement, with 
the exception of the two right-most peaks. These states 
lie in a region of the plane of winding ratios where not 
many POs are known and the action variables can be 
calculated only to correspondingly lower accuracy. 



VII. CONCLUSION 

The hydrogen atom in crossed electric and mag- 
netic fields presents long-standing challenges to both the 
dynamical-systems and the atomic-physics communities. 
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Figure 14: (a) The EBK spectrum as obtained from the quantization of the T£ (upper half, red) at E — —1.4 and F — 0.5 
compared to the exact quantum-mechanical spectrum (lower half, green). The first three n2-manifolds are labeled with (m, 
ri2, n-j). (b) Spectra at E = — 1.5, F = 0.5. 
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Figure 15: The quantization of the T3 1 family (blue) gives the 
two left-most peaks in the n% = 10 and 711 + ^3 = 10 (n = 10, 
q — 1) subseries showing a good agreement with the exact 
quantum spectrum (green). The other peaks result from the 
quantization of the T$ (red). The slight discrepancies of the 
peaks on the right are due to numerical inaccuracies in the 
interpolation procedure. 

In this work we have addressed both issues. From the 
viewpoint of dynamical systems we have demonstrated 
the power of periodic orbits when they are used as a 
probe to the intricacies of the geometrical and dynam- 
ical structures in phase space. Being the remnants of 
broken tori, periodic orbits can be used to establish a 
complete hierarchy of these phase space structures. At 
energies slightly below and slightly above the ionization 
saddle point of the crossed-fields hydrogen atom we es- 
tablished this hierarchy as follows. Three fundamental 
periodic orbits |4( that do not arise from the breakup of 
a higher-dimensional torus represent the 1-tori in the hi- 
erarchy Two of them serve as organizing centers for two 
families Tj 5 and TJj 1 of 2-torus POs. The periodic orbits 
located in the x-y symmetry plane form the T% family. 
The 2-torus POs themselves were found to be limiting 
cases of two families and T3 1 of 3-torus POs. 

Having established the hierarchy of POs, we calculated 
the individual action variables for the different families 



of POs. This knowledge provides the basis for the semi- 
classical calculation of the atomic spectrum using EBK 
quantization. Our results are in good agreement with the 
exact quantum mechanical spectrum. 

Because the classification of POs by winding num- 
bers relies only on the existence of a hierarchy of broken 
tori, which is a common feature in many non-integrable 
Hamiltonian systems, it will be applicable to other chal- 
lenging systems. In particular, previous experience on 
the hydrogen atom in a magnetic field |49| has shown 
that a semiclassical approximation that is derived in a 
near-integrable setting can successfully describe quantum 
states even deep in the mixed regular-chaotic regime. It 
is therefore to be expected that our quantization scheme 
will still be useful at appreciably higher field strengths 
than were considered here. 
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Appendix A: THE PERIODIC-ORBIT SEARCH 

Due to the Coulomb singularity in the Hamiltonian Q , 
a straightforward numerical integration of the trajecto- 
ries is unfeasible. To overcome this difficulty, we use 
the Kustaanheimo-Stiefel regularization jS^, ElJ an d in- 
tegrate in four- dimensional position and momentum co- 
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ordinates u and p. The equations of motion in KS co- 
ordinates with respect to a pseudotime parameter r are 
free of singularities (see, e.g., |47j). 

We describe the search for periodic orbits as a root- 
finding problem. A PO is identified by a starting point 
Pi = (u, p) i in phase space and by its pseudotime pe- 
riod r. Given an initial guess for Pj and t, we calculate 
the final point Pf of the trajectory by integrating the 
equations of motion for a time r and then use a Newton- 
like method to modify Pi (within the energy shell) and 
t so that Pf = Pi. Specifically, we use a Powell- hybrid 
root-finding method 0] for this task. 

The initial guesses required by the root-finder are ob- 
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